Cross-validation of active and passive microwave snowfall products over the continental United States

Surface snowfall rate estimates from the Global Precipitation Measurement (GPM) mission’sCoreObservatorysensors and theCloudSatradar are compared to those from the Multi-Radar Multi-Sensor (MRMS) radarcomposite product over the continental United States during the period from November 2014 to September 2020. Theanalysis includes the Dual-Frequency Precipitation Radar (DPR) retrieval and its single-frequency counterparts, the GPMCombined Radar Radiometer Algorithm (CORRA), theCloudSatSnow Profile product (2C-SNOW-PROFILE), and twopassive microwave retrievals, i.e., the Goddard Profiling algorithm (GPROF) and the Snow Retrieval Algorithm for GMI(SLALOM). The 2C-SNOW retrieval has the highest Heidke skill score (HSS) for detecting snowfall among the productsanalyzed. SLALOM ranks second; it outperforms GPROF and the other GPM algorithms, all detecting only 30% of thesnow events. Since SLALOM is trained with 2C-SNOW, it suggests that the optimal use of the information content in theGMI observations critically depends on the precipitation training dataset. All the retrievals underestimate snowfall ratesby a factor of 2 compared to MRMS. Large discrepancies (RMSE of 0.7–1.5 mm h21) between spaceborne and ground-based snowfall rate estimates are attributed to the complexity of the ice scattering properties and to the limitations of theremote sensing systems: the DPR instrument has low sensitivity, while the radiometric measurements are affected by theconfounding effects of the background surface emissivity and of the emission of supercooled liquid droplet layers

CloudSat-based assessment of GPM Microwave Imager snowfall observation capabilities

The sensitivity of Global Precipitation Measurement (GPM) Microwave Imager (GMI) high-frequency channels to snowfall at higher latitudes (around 60◦N/S) is investigated using coincident CloudSat observations. The 166 GHz channel is highlighted throughout the study due to its ice scattering sensitivity and polarization information. The analysis of three case studies evidences the important combined role of total precipitable water (TPW), supercooled cloud water,and background surface composition on the brightness temperature (TB) behavior for different snow-producing clouds. A regression tree statistical analysis applied to the entire GMI-CloudSat snowfall dataset indicates which variables influence the 166 GHz polarization difference (166∆TB)and its relation to snowfall. Critical thresholds of various parameters (sea ice concentration (SIC), TPW, ice water path (IWP)) are established for optimal snowfall detection capabilities. The 166∆TB can identify snowfall events over land and sea when critical thresholds are exceeded (TPW \u3e 3.6 kg·m−2, IWP \u3e 0.24 kg·m−2 over land, and SIC \u3e 57%, TPW \u3e 5.1 kg·m−2 over sea). The complex combined 166∆TB-TB relationship at higher latitudes and the impact of supercooled water vertical distribution are also investigated. The findings presented in this study can be exploited to improve passive microwave snowfall detection algorithms

SLALOM: An all-surface snow water path retrieval algorithm for the GPM microwave imager

This paper describes a new algorithm that is able to detect snowfall and retrieve the associated snow water path (SWP), for any surface type, using the Global Precipitation Measurement (GPM) Microwave Imager (GMI). The algorithm is tuned and evaluated against coincident observations of the Cloud Profiling Radar (CPR) onboard CloudSat. It is composed of three modules for (i) snowfall detection, (ii) supercooled droplet detection and (iii) SWP retrieval. This algorithm takes into account environmental conditions to retrieve SWP and does not rely on any surface classification scheme. The snowfall detection module is able to detect 83% of snowfall events including light SWP (down to 1 × 10−3 kg·m−2) with a false alarm ratio of 0.12. The supercooled detection module detects 97% of events, with a false alarm ratio of 0.05. The SWP estimates show a relative bias of −11%, a correlation of 0.84 and a root mean square error of 0.04 kg·m−2. Several applications of the algorithm are highlighted: Three case studies of snowfall events are investigated, and a 2-year high resolution 70°S–70°N snowfall occurrence distribution is presented. These results illustrate the high potential of this algorithm for snowfall detection and SWP retrieval using GMI

What Can We Learn from the CloudSat Radiometric Mode Observations of Snowfall over the Ice-Free Ocean?

The quantification of global snowfall by the current observing system remains challenging, with the CloudSat 94 GHz Cloud Profiling Radar (CPR) providing the current state-of-the-art snow climatology, especially at high latitudes. This work explores the potential of the novel Level-2 CloudSat 94 GHz Brightness Temperature Product (2B-TB94), developed in recent years by processing the noise floor data contained in the 1B-CPR product; the focus of the study is on the characterization of snow systems over the ice-free ocean, which has well constrained emissivity and backscattering properties. When used in combination with the path integrated attenuation (PIA), the radiometric mode can provide crucial information on the presence/amount of supercooled layers and on the contribution of the ice to the total attenuation. Radiative transfer simulations show that the location of the supercooled layers and the snow density are important factors affecting the warming caused by supercooled emission and the cooling induced by ice scattering. Over the ice-free ocean, the inclusion of the 2B-TB94 observations to the standard CPR observables (reflectivity profile and PIA) is recommended, should more sophisticated attenuation corrections be implemented in the snow CloudSat product to mitigate its well-known underestimation at large snowfall rates. Similar approaches will also be applicable to the upcoming EarthCARE mission. The findings of this paper are relevant for the design of future missions targeting precipitation in the polar regions

Chapter 12 - Recent advances and challenges in satellite-based snowfall detection and estimation

5noThe use of satellites for snowfall continuous monitoring and quantification on a global scale is fundamental to quantify water resources, and for understanding feedback mechanisms and interconnections in hydrology and climate. Spaceborne active and passive microwave sensors are particularly tailored to provide consistent measurements of falling snow thanks to their ability to penetrate clouds. However, remote sensing of snowfall remains among the most challenging tasks in precipitation science. Indeed, current state-of-the-art satellite products show large discrepancies in snowfall climatologies. The goal of this chapter is to provide an overview of the most recent studies dedicated to satellite-based snowfall detection and quantification by highlighting, in particular: (1) potentials and limitations of the most advanced spaceborne microwave sensors used for snowfall observation; (2) the main issues and challenges encountered in satellite-based retrieval of falling snow; and (3) recent advancements made in the scientific community in the attempt to overcome some of these issues. Concluding remarks include recommendations on future observation and retrieval strategies to improve detection and quantification of falling snow, with particular focus on higher latitudes.reservedmixedPanegrossi, Giulia; Casella, Daniele; Sanò, Paolo; Camplani, Andrea; Battaglia, AlessandroPanegrossi, Giulia; Casella, Daniele; Sanò, Paolo; Camplani, Andrea; Battaglia, Alessandr

The Passive Microwave Neural Network Precipitation Retrieval Algorithm for Climate Applications (PNPR-CLIM): Design and Verification

This paper describes the Passive microwave Neural network Precipitation Retrieval algorithm for climate applications (PNPR-CLIM), developed with funding from the Copernicus Climate Change Service (C3S), implemented by ECMWF on behalf of the European Union. The algorithm has been designed and developed to exploit the two cross-track scanning microwave radiometers, AMSU-B and MHS, towards the creation of a long-term (2000–2017) global precipitation climate data record (CDR) for the ECMWF Climate Data Store (CDS). The algorithm has been trained on an observational dataset built from one year of MHS and GPM-CO Dual-frequency Precipitation Radar (DPR) coincident observations. The dataset includes the Fundamental Climate Data Record (FCDR) of AMSU-B and MHS brightness temperatures, provided by the Fidelity and Uncertainty in Climate data records from Earth Observation (FIDUCEO) project, and the DPR-based surface precipitation rate estimates used as reference. The combined use of high quality, calibrated and harmonized long-term input data (provided by the FIDUCEO microwave brightness temperature Fundamental Climate Data Record) with the exploitation of the potential of neural networks (ability to learn and generalize) has made it possible to limit the use of ancillary model-derived environmental variables, thus reducing the model uncertainties’ influence on the PNPR-CLIM, which could compromise the accuracy of the estimates. The PNPR-CLIM estimated precipitation distribution is in good agreement with independent DPR-based estimates. A multiscale assessment of the algorithm’s performance is presented against high quality regional ground-based radar products and global precipitation datasets. The regional and global three-year (2015–2017) verification analysis shows that, despite the simplicity of the algorithm in terms of input variables and processing performance, the quality of PNPR-CLIM outperforms NASA GPROF in terms of rainfall detection, while in terms of rainfall quantification they are comparable. The global analysis evidences weaknesses at higher latitudes and in the winter at mid latitudes, mainly linked to the poorer quality of the precipitation retrieval in cold/dry conditions

A Machine Learning Snowfall Retrieval Algorithm for ATMS

This article describes the development of a machine learning (ML)-based algorithm for snowfall retrieval (Snow retrievaL ALgorithm fOr gpM–Cross Track, SLALOM-CT), exploiting ATMS radiometer measurements and using the CloudSat CPR snowfall products as references. During a preliminary analysis, different ML techniques (tree-based algorithms, shallow and convolutional neural networks—NNs) were intercompared. A large dataset (three years) of coincident observations from CPR and ATMS was used for training and testing the different techniques. The SLALOM-CT algorithm is based on four independent modules for the detection of snowfall and supercooled droplets, and for the estimation of snow water path and snowfall rate. Each module was designed by choosing the best-performing ML approach through model selection and optimization. While a convolutional NN was the most accurate for the snowfall detection module, a shallow NN was selected for all other modules. SLALOM-CT showed a high degree of consistency with CPR. Moreover, the results were almost independent of the background surface categorization and the observation angle. The reliability of the SLALOM-CT estimates was also highlighted by the good results obtained from a direct comparison with a reference algorithm (GPROF)

Rainfall microphysical characterization over the Mediterranean area during the GPM era

The NASA/JAXA Global Precipitation Measurement (GPM) Core Observatory (CO) was launched on 27th February, 2014. It carries, for the first time, a Dual-frequency Precipitation Radar (DPR) designed to provide insights into the 3-D structure of precipitating clouds and rain intensity by using its Ka- and Ku-band frequencies. In addition to characterize the 3-D structure of precipitation, the DPR is used as calibrator for the GPM Microwave Imager (GMI). Single-frequency (SF) (both Ku- and Ka-only) and double-frequency (DF) based products provide, among the others, particle-size distribution (PSD) parameters (namely the mass-weighted mean diameter Dm and the normalized intercept parameter Nw), as well as precipitation rates. This book chapter focuses on reliability of the PSD parameters (limited to rainfall events, in this case we'll talk about of DSD - drop size distribution) over the Mediterranean area by taking as reference the DSD parameters estimated by ground-based radar measurement. Before of this, an inter-comparison between the SF and DF DPR outputs considering five years of data from GPM-CO mission is carried out. The goal is to investigate the reliability of SF-based products by assessing their quality compared to the DF-based ones, treated as a reference

A 4-Year Climatological Analysis Based on GPM Observations of Deep Convective Events in the Mediterranean Region

Since early March 2014, the NASA/JAXA Global Precipitation Measurement Core- Observatory (GPM-CO) satellite has allowed analysis of precipitation systems around the globe, thanks to the capabilities of the GPM Microwave Imager (GMI) and Dual-Frequency Precipitation Radar (DPR). In this work, we demonstrate how GPM-CO measurements obtained from 4 years of observations over the Mediterranean area can be used as an extremely effective tool to study the main climatological characteristics of the most intense Mediterranean storm structures. DPR and GMI-based Precipitation Features (PFs) parameters are used as proxies of the vertical structure and microphysical properties of these events, and their statistical distribution is analyzed to identify extremes. The analysis of annual, seasonal and geographical distribution of the identified deep convective systems highlights substantial differences in their diurnal cycle and in the distribution between land-sea and summer-winter. There is a general shift of the convective systems from the south (mostly over the sea) in the cold season, to the north (mostly over land) in the warm season. The analysis shows also that the inferred convective intensity is not always related to heavy precipitation. Known DPR and GMI-based criteria were adopted to identify overshooting top events and potential hailstorms, identify extreme deep convection signatures, like those observed for tropical and subtropical systems, and the most intense occur mostly over the sea. Although the analysis is limited to four years, the results show that the GPM-CO offers unprecedented measurements to identify and characterize extreme weather events in the Mediterranean region, with unique potentials for future long-term climatology and interannual variability analysis

Evaluation of the Sensitivity of Medicane Ianos to Model Microphysics and Initial Conditions Using Satellite Measurements

Tropical-like cyclone (TLC or medicane) Ianos formed during mid-September 2020 over the Southern Mediterranean Sea, and, during its mature stage on days 17–18, it affected southern Italy and especially Greece and its Ionian islands, where it brought widespread disruption due to torrential rainfall, severe wind gusts, and landslides, causing casualties. This study performs a sensitivity analysis of the mature phase of TLC Ianos with the WRF model to different microphysics parameterization schemes and initial and boundary condition (IBC) datasets. Satellite measurements from the Global Precipitation Measurement Mission-Core Observatory (GPM-CO) dual-frequency precipitation radar (DPR) and the Advanced Scatterometer (ASCAT) sea-surface wind field were used to verify the WRF model forecast quality. Results show that the model is most sensitive to the nature of the IBC dataset (spatial resolution and other dynamical and physical differences), which better defines the primary mesoscale features of Ianos (low-level vortex, eyewall, and main rainband structure) when using those at higher resolution (~25 km versus ~50 km) independently of the microphysics scheme, but with the downside of producing too much convection and excessively low minimum surface pressures. On the other hand, no significant differences emerged among their respective trajectories. All experiments overestimated the vertical extension of the main rainbands and display a tendency to shift the system to the west/northwest of the actual position. Especially among the experiments with the higher-resolution IBCs, the more complex WRF microphysics schemes (Thompson and Morrison) tended to outperform the others in terms of rain rate forecast and most of the other variables examined. Furthermore, WSM6 showed a good performance while WDM6 was generally the least accurate. Lastly, the calculation of the cyclone phase space diagram confirmed that all simulations triggered a warm-core storm, and all but one also exhibited axisymmetry at some point of the studied lifecycle